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ABSTRACT 

In the course of the follow-up multiwavelength study of a deep radio survey 
we have discovered that the milli-Jansky radio source PDFJ011423 is a low- 
redshift {z = 0.65) extremely red galaxy (ERG) with K = 15.3, R — K = 5.8 and 
J — K = 3A. Optical, infrared and radio photometry, together with optical and 
near-infrared spectroscopy, reveal a heavily obscured galaxy (^4^=5-6, from the 
observed Balmer decrement) undergoing vigorous star formation and presenting 
an active galactic nucleus (AGN). PDFJ011423 is a representative member of the 
dusty ERG population, providing a local counterpart for studying more distant 
ERGs. 

Subject headings: galaxies: individual (PDFJ011423) — galaxies: photometry — 
galaxies: active — galaxies: starburst — infrared: galaxies — radio continuum: 
galaxies 
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1. INTRODUCTION 

Deep multi-wavelength near-infrared and optical surveys over the last decade have led 
to the discovery of a population of galaxies with very red optical-infrared colors ("extremely 
red galaxies" , hereafter ERGs; e.g. Elston, Rieke & Rieke 1988; McCarthy, Persson & West 
1992; Hu & Ridgway 1994; Thompson et al. 1999; Daddi et al. 2000). ERGs, defined as 
galaxies with R — K > appear to fall into two sub-classes: (1) high redshift {z > 1) 
evolved eUiptical galaxies with intrinsically red spectral energy distributions (SEDs) (e.g. 
Spinrad et al. 1997; Stiavelli et al. 1999; Soifer et al. 1999) and (2) highly obscured dusty 
galaxies undergoing starburst activity (e.g. Gimatti et al. 1998; Dey et al. 1999; Small et 
al. 1999). The predominance in the ERG population of the latter would favor hierarchical 
galaxy formation scenarios (White & Frenk 1991), while the former would support early 
formation of galaxies through monolithic collapse (Larson 1975). Therefore, the study of 
ERGs can help constraining existing models for formation and evolution of galaxies. 

The observed very red colors of those ERGs which are classified as elliptical galaxies are 
attributed predominantly to the large K-correction arising from their high redshift. On the 
other hand, star-forming galaxies will owe their redness predominantly to dust extinction. 
This view implies that extremely red ellipticals will occur only at high redshift {z > 1) while 
extremely red galaxies at 2; < 1 will exhibit unusually heavy obscuration of either stellar or 
AGN emission, or both. This is expected to be the reason for the small space density found 
for ERGs in the local Universe. 

This letter reports the discovery of an ERG, PDFJ011423, at z ^ 0.65. Multi-waveband 
photometric and spectroscopic observations reveal evidence of both star formation and AGN 
activity in this galaxy, which can be regarded as a "local" template for the study of this class 
of ERG at higher redshifts. Throughout this paper we adopt Ho — 65 /i65 kms~^ Mpc~^ and 
go = 0.5. 

2. OBSERVATIONS 

PDFJ011423 {a = 0l'^14"*23^ 5 = -45°34'30", J2000) was first noted as a faint 1.4 GHz 
radio source (5'i.4ghz = 1.67 mJy) in the Phoenix Deep Survey (Hopkins et al. 1998, 1999). 

Aperture photometry of the optical counterpart of the source gives i?=21.1 and V=22.7 
(Georgakakis et al. 1999), locating PDFJ011423 close to the magnitude limit of the follow-up 
Two Degree Field multi- fibre spectroscopic survey (Georgakakis et al. 1999). Nevertheless, 
a strong emission line, identified as [OII]3727 at a redshift z = 0.65, was detected. 

In the light of the known strong correlation between the far-infrared (FIR) and 1.4 GHz 
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luminosities for star-forming galaxies, these clues were sufficient to schedule PDF JOl 1423 

for observation by the Infrared Space Observatory (ISO), using ISOCAM (7 and 15 /xm) 
and ISOPHOT (90 /xm). The /S'O observations (J. Afonso et al, in preparation) revealed 
the presence of a relatively bright source at 7 and 15 ;um, with fluxes of 4.1 and 7.6 mJy, 
respectively. The object was also detected by ISO at 90 fim (4(7 level), with a flux of 260 mJy. 
Subsequent ncar-lR photometry, using the New Technology Telescope (NTT) and the CTIO 
1.5m telescopes, showed this galaxy to have K—15.S and J=18.4. PDFJ011423 is thus 
classified as an ERG with R — K — 5.S and J — K — 3.1. 

The closeness and the relatively bright optical luminosity of PDFJ011423 allows detailed 
spectroscopy of this object at both optical and ncar-lR wavelengths. The optical spectro- 
scopic observations were obtained at the Very Large Telescope (VLT) using the FORSl 
instrument in longslit mode. Three dithered observations of 15 minutes each were made 
using a 1.6" wide slit and a low resolution grism (150I-I-17), covering the wavelength range 
AAO. 6-1.1 //m. Bias subtraction, flatfielding and the combination of the dithered observa- 
tions were performed before extracting the final spectrum using the software package IRAF. 
The wavelength scale was defined by fitting a third order polynomial to the lines of a He-Ar 
calibration spectrum. The spectrum of the standard star GD50 was taken for the fiux cali- 
bration, which agrees with the optical photometry of PDFJ011423 within 15%. Near-infrared 
spectroscopy was obtained using SOFT on the ESO New Technology Telescope. Twelve ob- 
servations of 270 seconds each, dithered along a 1.0" sht, were obtained. The low resolution 
grism covering the range AAO. 95-1. 64 //m was used. The dithered observations were flat- 
fielded and combined before the spectrum was extracted. For the wavelength calibration a 
Xe lamp spectrum was used. The photometric calibration was performed using the J-band 
magnitude of the galaxy. The flnal calibrated spectra are presented in Figure 1 with the 
respective line measurements, based on Gaussian flts to the emission lines, given in Table 1. 



3. NATURE OF PDFJ011423 

We approach the interpretation of the observed spectral energy distribution (SED) of 
PDFJ011423 (Figure 2) by attempting flrst to model its continuum emission, and then 
checking for consistency with the observed spectral line emission. The SED of a normal 
elliptical galaxy can only reproduce the observed very red colors of PDFJ011423 at z > 1. 
This contradiction with the known redshift, combined with the presence of strong emission 
lines, eliminates PDF JOl 1423 as an elliptical galaxy. 

The high luminosity of the galaxy at rest-frame 60 fim suggests the presence of a dusty 
star burst, while not excluding the possibility that some of the near-IR and mid-IR luminos- 
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ity arises from nuclear activity. Accordingly, we have matched the observed fluxes by com- 
bining the model SEDs of a starburst occurring in evolving, dusty giant molecular clouds 
(GMCs) (Efstathiou, Rowan- Robinson & Siebenmorgen 2000) and of a dust-enshrouded 
AGN (Rowan- Robinson 1995). Although it is possible that the two energising phenomena 
are physically linked (e.g. Rowan- Robinson 1995), here we simply superimpose the SEDs 
predicted by the two models. 

The model fit, presented in Figure 2, results from a minimization to the observed 
SED, varying within the domain of the parameters of the original starburst and AGN models. 
The more powerful component, responsible for 76% of the infrared luminosity (ivi-iooo^jm), is 
a model of a dusty starburst of age 57Myr, with optical opacities of the GMCs in the range 
Ty = 3 — 200. The properties of this component are determined primarily by the constraints 
imposed by the FIR and optical emission. 

The starburst component alone is, however, unable to fit the mid-IR data. Therefore, 
an extra contribution in the form of emission from hot dust (T ~ 160 — 1600 K) heated by a 
putative AGN is required. This component represents 24% of the total infrared luminosity. 
Acceptable fits could also be obtained using other AGN models, for example by allowing for 
emission from a dusty torus surrounding the central (black hole) source (Efstathiou, Hough 
& Young 1995). Both classes of models lead to similar descriptions of the object in the 
context of this paper. 

The best fitted model SED in Figure 2 predicts an infrared luminosity of Li_iooo/im = 
7.1x10^^ h^i Lq, which classifies PDFJ011423 as an ultra-luminous infrared galaxy (ULIRG) . 
The possibility that the high luminosity is due to gravitational lensing is unlikely, given the 
relatively low redshift of the object. 

The existence of only one data point in the FIR/sub-mm part of the SED makes an 
estimate of the dust mass for PDF JOl 1423 highly uncertain. However, the fitted SED in 
this wavelength range is compatible with an optically thin thermal dust emission spectrum 
(emissivity index of 1.5) with a temperature of Tjust = 33 ± 3K, which corresponds to a 
dust mass of Mdust ~ 8 - 30 X 10*^ h^i Mq (Hildebrand 1983). Consistent with the extreme 
reddening of the galaxy, this value is an order of magnitude higher than the dust mass of 
10^ - 10^ Mq found for local ULIRGs (Sanders & Mirabel 1996) and comparable to the 
higher values found in the sample of PG quasars of Haas et al. (2000), for example. 

The two models used to fit the observed SED do not predict the radio continuum 
luminosity or spectral index. To interpret the radio data, we appeal to the remarkably 
tight empirical correlation between the FIR and radio continuum luminosity observed in 
star- forming galaxies (Helou, Soifer & Rowan- Robinson 1985; Condon 1972). The favored 
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explanation of the correlation holds that the same massive stars warm the FIR-emitting 
dust and energise, through supernova explosions, the relativistic electrons responsible for 
the radio continuum. Given the FIR luminosity of PDFJO 11423, this correlation predicts a 
1.4 GHz flux in the range 1.3 - 3.3 mJy compared with the measured value of 1.67 mJy. This 
agreement supports the possibility that most of the radio emission comes from star formation 
processes. It should however be noted that a similar radio/FIR correlation holds for samples 
of radio-quiet quasars, perhaps as a consequence of links between star formation and black 
hole feeding rates. The relatively flat radio spectral index observed here, aff — 0.16 (S,^ oc 
i/"""), while being frequently hnked to quasar emission (e.g. Webster et al. 1995), has also 
been observed in compact starburst nuclei of ULIRGs (Sopp & Alexander 1991; Crawford et 
al. 1996). The available radio images do not have sufficient spatial resolution to reveal any 
structure in PDFJ011423 and cannot aid to clarify the origin of the radio emission. 

We now consider the spectral properties of PDFJ011423, using optical and near-IR 
line diagnostics from Figure 1. Since the Hq;-|-[NII] blend is not resolved by the NIR spec- 
trum, we assume the Ha/ [Nil] ratio to have a value between 1 (typical of Seyfert galaxies, 
e.g. Veilleux & Osterbrock 1987) and 2 (corresponding to star-forming galaxies, Kennicutt 
1992) in order to separate the two contributions. With the de-blended H« line intensity, the 
observed Balmer line ratio becomes Ha/H/? = 15 — 20. Assuming Case B Balmcr recombina- 
tion (adopting the intrinsic value for Ha/H/S as 2.85 for star- forming galaxies (Brocklehurst 
1971) and 3.1 for narrow line AGNs (Veilleux & Osterbrock 1987)) and a standard Galactic 
extinction curve (GardeUi, Glayton & Mathis 1989, using Ry — 3.1), we estimate an op- 
tical dust extinction of Ay — 5 — 6 mag. This large extinction is hkely to be the source 
of the extreme red color in PDFJ011423. The absence of the 4000 A break in the opti- 
cal spectrum, indicating an optical continuum dominated by very young stars, is consistent 
with the inference that this is a vigorously star-forming galaxy. The measured optical line 
widths are not significantly larger than the instrumental resolution, placing an upper limit 
of approximately 500 km s^^ on their intrinsic value. In the NIR spectrum, the Hq;+[NII] 
blend presents a width value of (FWHM) ~ 2000 km s~^ (after correction for instrumental 
resolution) . Given the low resolution of the spectrum it is not clear if this shows an intrinsic 
broad line component for Ha or is due to significant contribution from [Nil]. 

The diagnostic fine ratios (Rola, Terlevich & Terlevich 1997; Veilleux & Osterbrock 
1987) for PDF JOl 1423, after correction for dust obscuration, are consistent with the two 
component model. While ratios involving the [OII]3727 fine ([OII]3727/H/3 = 8.0), indicate 
excitation by an AGN, the other line ratios ([OIII]5007/H/3=2.9, [SII]/Hq;=0.2-0.3) and the 
non-detection of [01] 6300 are indicative of a starburst. This kind of ambiguous spectral line 
classification has been linked previously to galaxies composite in nature, i.e., hosting both a 
starburst and an AGN (HiU et al 1999). 
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The presence of an AGN in addition to the starburst component complicates the estima- 
tion of the star formation rate (SFR). The most rehable measure is from its 60 /xm luminosity 
(2.4 X 10^^ /igg^WHz"^) where the model implies that the dominant contribution is due to 
the starburst. Using the appropriate calibration (Cram et al. 1998) for an initial mass func- 
tion (IMF) of the form ^(M) ~ M^^'^ with stellar masses between 0.1 and 100 M©, we 
estimate a SFR for massive stars of SFReo^m (M > 5 M0)= 470 h^^ M© yr'^ Alternatively, 
from the observed radio power (Pi.4ghz = 1-5 x 10^^ /igg^WHz"^) we deduce SFRi.4ghz 
(M > ^Mq)— 377 h^^ Mq yr~^, consistent with the value estimated from the 60 jjim lumi- 
nosity. On the other hand, the extinction-corrected Ha luminosity (Lhq = 0.9 — 2. 6x1 0^^ K^^ 
W) implies SFRHa(M > 5 M0)= 620 - 1730 /igg^ M© yr^S the large range being due to the 
unresolved contribution of [Nil] to the observed IIa+ [Nil] blend. This excess over the SFR 
estimated from FIR and radio luminosities may be due to the AGN contribution. 

Three other dusty ERGs have been spectroscopically confirmed and are relatively well 
studied: HRIO, displaying strong starburst activity at 2; = 1.44 (Dey et al. 1999); ISOJ1324- 
2016, a. z = 1.50 galaxy which hosts a dusty quasar (Pierre et al. 2001) and EROJ164023, 
which dX z = 1.05 shows star-forming activity with a possible weak AGN (Smith et al. 2001). 
A brief comparison of some of the observed quantities in PDFJO 11423 and these three other 
dusty ERGs is listed in Table 2. Although the non detection of star formation in ISOJ1324- 
2016 can be due to the lack of FIR/sub-mm observations, it seems clear that different degrees 
of AGN and starburst activity do exist in the dusty ERG population. Given the lower redshift 
of PDFJ011423, this source offers the best opportunity to study the interplay between the 
two phenomena in the dusty environments of these galaxies. 

The ratio of the radio power to the dust-corrected optical luminosity of PDFJ011423, 
Pi.a/Lr ~ 1, confirms that PDFJ011423 is radio quiet. It seems hkely that this source is 
linked to the radio-quiet counterparts of the red quasar population discovered by Webster 
et al. (1995) and also to the ERGs recently discovered among the population responsible for 
the hard X-ray background (Hasinger et al. 2001). Future work will show if the intense star 
formation activity present in PDF JOl 1423 is also a common feature to the above mentioned 
obscured AGN populations. 

All these different studies suggest the existence of an important population of dusty 
ERGs, powered by heavily obscured starbursts and/or AGNs, which is now starting to be 
observed. Given their extreme nature, the study of the dusty ERGs will hold fundamental 
clues to the understanding of galaxy evolution, by revealing valuable information on the 
hidden star formation and AGN activity in the universe. 

We are grateful to A. Verma for adapting and running her SED fitting algorithm for this 
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Table 1. Line measurements for PDF JOl 1423 



Line 


Aobs 




Fobs 


W^A,rest 




(A) 




(10-^6 erg s-^ cm-2) 


(A) 


[OII]AA3726,3729 


6159 ± 


1 


8.0 ± 0.3 


90 ± 6 


H/3 


8033 ± 


4 


1.0 ± 0.2 


10 ± 4 


[OIII]A4959 


8195 ± 


4 


2.0 ± 0.8 


19 ± 11 


[OIII]A5007 


8274 ± 


4 


3.7 ± 0.7 


32 ± 10 


Hq;+[NII]'^ 


10858 ± 


11 


28.7 ± 2.5 


217 ± 47 



^Measurements from a single Gaussian fit to the observed blend. 
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Table 2. Comparison between the observed properties of dusty ERGs 



PDFJ011423 



HRIO^ ISOJ1324-20162 EROJ164023 



z 

K mag 
R-K 
I -K 
J-K 
Lna (W) 

Lna.corr (W)^ 
J-'IO /iin— 2 cm (Lo) 

SFR (M> 5M0, Meyr-i) 
Mdust (Mo) 
Av (mag) 



0.65 
15.3 
5.8 

3.1 

2.0 - 2.6 X 10^5 
0.9 - 2.6 X 10^7 
6.5 X 10^2 
~ 470 
8 - 30 X 10« 
5-6 



1.44 
18.4 

7.8 
2.6 
2.6 X 10^5 

4.0 X 10^2 
200 - 400 
~ 4 X 10^ 
2-4.5 



1.50 
17.5 
> 5.2 
4.9 



1.05 
17.6 
5.9 
4.5 
2.1 

0.7 X 10^5 

< 2 X 10^2 
~ 6 - 400 



4-7 



^Ha luminosity corrected for dust extinction. 



Note. — All the values are given for the cosmology adopted in this paper. 



References. — (1) Dey et al. 1999; (2) Pierre et al. 2001; (3) Smith et al. 2001 
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Fig. 1. — Optical and near-infrared spectra of PDFJ011423 at z = 0.65. The optical spec- 
trum was obtained at the Very Large Telescope (VLT) using FORSl while the near-infrared 
spectrum was obtained at the New Technology Telescope (NTT) using SOFI. 
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log (i/[Hz]) 

Fig. 2. — Observed SED of PDF JOl 1423 compared with the fitted starburst (long-dashed 
fine) and AGN dust emission (short-dashed fine) models. Crosses and stars denote, re- 
spectively, the optical (V and i?-bands) and near-infrared (J and A'-bands) measurements. 
Detections at mid and far-infrared are from ISO (filled diamonds and square) while upper 
limits (3(7, open circles) are shown for co-added observations from the Infrared Astronomical 
Satellite (IRAS) at Aobs = 100, 60, 25, and 12 //m, estimated using the SCANPI procedure 
at the Infrared Processing and Analysis Center. Filled circles are the radio detections at 
1.4 and 2.4 GHz. The shaded region represents the prediction of the radio fiux based on 
FIR/radio correlation described in the text. 



